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ABSTRACT 

The observation of the polarization emerging from a rotating star at different 
phases opens up the possibility to map the magnetic field in the stellar surface 
thanks to the well-known Zeeman Doppler Imaging. When the magnetic field is 
sufficiently weak, the circular and linear polarization profiles locally in each point 
of the star are proportional to the first and second derivatives of the unperturbed 
intensity profile, respectively. We show that the weak-field approximation (for 
weak lines in the case of linear polarization) can be generalized to the case of a 
rotating star including the Doppler effect and taking into account the integration 
on the stellar surface. The Stokes profiles are written as a linear combination 
of wavelength-dependent terms expressed as series expansions in terms of Her- 
mite polynomials. These terms contain the surface integrated magnetic field and 
velocity components. The direct numerical evaluation of these quantities is lim- 
ited to rotation velocities not larger than 8 times the Doppler width of the local 
absorption profiles. Additionally, we demonstrate that, in a rotating star, the cir- 
cular polarization flux depends on the derivative of the intensity flux with respect 
to the wavelength and also on the profile itself. Likewise, the linear polarization 
depends on the profile and on its first and second derivative with respect to the 
wavelength. We particularize the general expressions to a rotating dipole. 

Subject headings: Stars: magnetic fields, rotation — Techniques: polarimetric - 
Methods: analytical 
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Introduction 



The present knowledge of solar and stellar magnetism derives essentially from the 
interpretation of polarization in spectral lines due to the Zeeman effect^. In the presence 
of a weak magnetic field, the radiative transfer equation for polarized light can be solved 
analytically and it is known as the "weak field approximation". Although very simple, 
it relies on a set of strong assumptions that we discuss later. However, the weak field 
approximation is applicable in many scenarios and gives very good results for the inference of 
magnetic fields as compared with more elaborate methods. In solar physics, it is at t he heart 



of the success of a 



1995 



Spirock et al 



arge 



number of synoptic magnetographs, like those of Big Bear (Varsi 



200ll ). It has been also used recently to produ ce vector magnetogram s 



with the Imaging Magnetograph Ex periment (IMaX) in strument (IMartmez Pillet et al. 



201 ll ) onboard the Sunrise balloon (ISolanki et al.ll2010l ). In stellar spectropolarimetry. 



the weak fie 



Donati et al. 



d app roximation is at the base of the least-squares deconvolution (LSD; 



19971 ). the most successful technique used to detect and measure magnetic 



fields in solar-type s tars, or in some r ecent works on c entral stars of planetary neb ulae 



( iJordan et al. 



2005 



pulsat ing stars ( 



Leone et al. 



Silvester et al. 



stars ( jWade et al. 



20001; 



2009). 



20111 ). white dwarfs 



rot su b dwarfs ( lO'Toole et al. 



Bagnulo et al. 



fjAznar Cuadrado et al. 



2004J) 



20021 ). 



20051 ) and Ap and Bp 



Most of the scientific cases commented so far deal with unresolved structures which 
means that the weak field approximation has to be under stood in terms of fluxes o f the 



Stokes parameters instead of specific intensities. Recently, 



Martinez Gonzalez et al. 



f l2012h 



derived the weak field expressions in terms of fluxes, particularizing to the case of the 



1 There are some notable exceptions as, for example, the study of the quiet Sun magnetism 



by means o 
effect (e.g. 



the interpretation o f scattering polarization and its modification by the Hanle 



Trujillo Bueno 



2011 



for a recent review) 
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stellar dipole. They neglected rotation to simplify the equations and to derive analytical 
expressions for the inference of the magnetic field. When rotation is taken into account (in 
other words, there is a correlation between the magnetic field and the line-of-sight velocity), 
the problem becom es much more challenging and it is usually solved numerically (e.g., 



Petit fc Wade 



20121 . for a recent effort). However, we show in this paper that it is possible 
to obtain analytical expressions for the Stokes vector fluxes in the weak field approximation. 
There are two interests on this effort. First, having an analytical expression for the Stokes 
flux might be of help to investigate the interplay between rotation and the magnetic field. 
Second, it can be used to speed up inversion codes based on the weak-field approximation 



because no numerical integration is neede d. This will fac i 



inference codes, like the one developed by 



itate the application of Bayesian 



Petit fc Wadd ( 120121 ). that are based on Markov 



Chain Monte Carlo methods that need to carry out the synthesis thousands of times to 
sample the posterior distribution. 

This is the first paper of a series dealing with analytical expressions for the fluxes of 
the Stokes parameters in the weak field limit. This paper deals with the expressions for a 
rotating star with a magnetic field on its surface. The results depend on the correlation 
between the magnetic field and the velocity and are general for any rotation profile or 
any magnetic field configuration. We particularize the equations to the case of a rotating 
magnetic stellar dipole. 



2. The weak field approximation for a star 

Let us consider a star with arbitrary velocity and magnetic fields on its surface. 
We define both fields for an arbitrary point r in the surface using the reference frame 
represented in Fig. [TJ The line-of-sight (LOS) is along the Z axis, while we choose the axis 
X as the one defining the reference direction for positive Stokes Q. 
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Fig. 1. — This figure shows the reference systems used to define the Stokes parameters in 
this work. We consider that the observer is along the Z axis. The direction X defines the 
reference for Q > 0. The vector r defines any position on the stellar surface. 

Our aim is to compute the local Stokes profiles emerging from every point of the visible 
surface and compute its integral. As explained in the introduction, a simple (although 
subject to approximations) way to relate the emergent Stokes profiles with the magnetic 
field vector is the weak-field approximation. This approximation (combined with the 
weak-line approximation for linear polarization) leads to an analytical solution to the 
radiative transfer equation. It is based on the following three fundamental assumptions: 
i) the magnetic field vector is constant along the line of sight on the region of the 
stellar atmosphere where the spectral line is formed, ii) the magnetic field intensity is 
sufficiently weak so that the Zeeman splitting is smaller than any other line broadening 
mechanism and iii) the LOS velocity and any broadening mechanism are constant with 
height in the line formation region. Under these simplifying conditions, the magnetic 
field can be considered as a perturbation to the zero-field case. This approximation is 
not strictly valid in general, but it is a good approximation whenever the magnetic field 
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strength is weak in the line formation region in the absence of strong gradients along 
the LOS. Quantit atively, the approximation holds whene ver gAXs <C AA^, where AA# 



(see Eq. (3.14) in 



Landi Pegl'Innocenti fc Landolfi 



20041 ) is the Zeeman splitting that is 



proportional to the field strength, A Ad is the dominant broadening mechanism and g 
represents the effective Lande factor of the line, which quantifies the magnetic sensitivity 
of the line and that only depends on the quantum numbers of the transition (see, e.g., 



Landi DeglTnnocenti fc Landolfi 



20041 ). From its definition, it is clear that the weak field 



regime occurs at different field strengths for different spectral lines, for different stellar 
objects and for different spectral resolutions. For instance, for a magnetically sensitive iron 
spectral line with g = 3 at 5000 A with a standard thermal broadening at T = 6500 K, 
the weak-field regime is applicable for fields below ~600 G if the spectral resolution of the 
spectrograph is smaller than the line broadening. 

Under the previous simplifying assumptions, and to first order in AA#, the intensity 
profile /(A) of a spectral line formed in a weak magnetic field is fully insensitive to the 
magnetic field. In other words, it fulfills the s tandard radiative transfer equation in the 



absence of a magnetic field (e.g., 



Mihalas 



19781 1 . Likewise, at first order in A As, the circular 



polarization profile, i.e., the Stokes V^(A) profile for a given spectral line, has the following 
expression: 



V(X) = -CXlgB, 



9_m 

dX 



where the symbol B z stands for the longitudinal component of the magnetic field (the 
projection along the LOS in the reference system shown in Fig. [T]). The symbol A stands 
for the wavelength while Ao is the central wavelength of the transition. The constant 
C = 4.67x lO-^G^A -1 . 

At first order in AA#, linear polarization is zero. In order to obtain an expression 
for the Stokes profiles characterizing linear polarization we have to expand the radiative 



-7- 



transfer equation to second order in AXb and assume that the spectral line is not saturated 
and weak. Under these assumptions, a good approximation to the intensity profile emerging 



from the stellar atmosphere is to consider it, in general, proportiona 



(see Eqs. (9.84) and (9.86) in 



Landi Degl'Innocenti fc Landolfi 



to a Voigt function 



20041 ) . For the sake of 



simplicity, we consider the case of a Gaussian line profile. All our results can be derived 
for a Voigt profile but the final expressions are more complicated without any fundamental 
difference. This assumption somehow reduces the scope of applicability of our results for 



linear polarization because the line shape self- similarit y is only valid for 



depth is smaller than ~40% according to the results of 



Kochukhov et al. 



ines whose line 



(120 10h . 



With this assumption, the weak-field approximation f or the Stokes Q(\) and U(X), de 



fined in the XYZ reference system sh own in Fig. [TJ are (see 



2004 



Martinez Gonzalez et al. 



20121 ): 



Landi Degl'Innocenti fc Landolfi 



QW 



^xWr 2 P 2 ^ 2 -f(A) 



U(X) = — —\qGB x B v ^ 2 - , 



(2) 



where the symbol G plays the role of the effective Lande factor for linear polarization 



and quantifies the sensitivity of linear polarization to t 



re magnetic field. Again, it is on ly 



a fun ction of the quantum numbers of the transition f lLandi Degl'Innocenti fc Landolfi 



20041 ) . The symbols B x and B y stand for the components of the magnetic field vector along 
the X and Y axes, respectively. The Stokes profiles in any other reference system that 
i s defined as a rotation around the Z a xis by an angle 7 are given by the usual rotation 



( Landi Degl'Innocenti fc Landol 



2004|): 



Q ncw (A) = cos(2 7 )Q(A) + sin(2 7 )f/(A) 
EWA) = -sin(2 7 )Q(A) + cos(2 7 )f/(A), 



(3) 



where the quantities with the subindex "new" refer in this case to the new reference frame. 
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When observing unresolved objects, the detected polarized radiation is obtained as an 
integration in the plane of the sky of the individual Stokes parameters at each point of the 
stellar surface (E), all of them referred to a common reference system. We describe the 



stellar disk using the polar coordin ates p and 8, where p is norma 



so that its range is from to 1 (see 



Martinez Gonzalez et al. 



ized to the stellar radius, 



20121 ). Therefore, the average 



of any arbitrary function h(p, 9) on the surface £ will be given by: 

p 1 r-l'K 

I dEh(p, 6)= pdp / d0h(p, 6) 
Jo Jo 




(4) 



As a consequence, it is obvious that the observed values of the Stokes flux vector depends 
on the surface distribution of the magnetic field, on the centre-to-limb variation (CLV) of 
the radiation, and on the Doppler effect due to the rotation of the star. 

Follow ing; standard approximat i ons a nd assuming that the star has azimuthal 



symmetry, 



Martinez Gonzalez et al 



( 120121 ) parameterized the intensity at each point of 



the star as a product of two functions: Iq(\), that only depends on the wavelength (the 
intensity at the disc centre since they neglect rotation), and f( u), the CLV t 
on the position on the disk and is wavelength independent (e.g. 



Claret 


2000; 


Cox 


2000) 



f(p) = (1 — a — b + ap + bp 2 ). 



(5) 



The CLV is given in terms of p = cos = y/l — p 2 , where is the astrocentric angle 
between the normal to a point in the stellar surface and the line of sight (see Fig. [T]). The 
parameters a and b have values between and 1 and are supposed to be constant along 
the spectral line (although they can vary from line to line). With this assumption, we also 
force the line depth to be //-independent. Note that the values of a and b have to fulfill 
the condition that the intensity has to be positive. This simple parameterization allowed 
them to express the weak field equations in terms of the observed fluxes, thus mimicking 
the expressions for the specific Stokes vector. This was possible because the intensity flux 
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is assumed to be the product of the wavelength-dependent intensity at disk center and 
the wavelength independent CLV. Using this relationship, the weak-field and weak-line 
approximation for the integrated polarized flux remains formally the same than for the 
specific Stokes vector but the components of the magnetic field appear weighted by the 
CLV law. 

When the velocity field is included, the situation turns out to be much more complicated 
because the intensity itself depends on the local projection of the velocity along the Z axis 
(v z ), which subsequently depends on the position along the stellar disk in a potentially 
complicated way. In this case, the fluxes of the Stokes vector are expressed as: 

Fj(A) = f pdp f dd I{\,v z {p,6),n) 
Jo Jo 

f q (x) = -^xt G J\d P f\e [Bl( P ,e)-Bl( P ,e)] *«^(p,0U) 

Fu(X) = -^XtG jT pdp jT d6 B x (p, 6)B y ( P , 6) ^ ^ & 
F V {X) = -CXlg f\dp J\9 B z (p,9) d J^^ , 



(6) 



The immediate consequence of the dependence on v z (p,6) is that the integrals cannot 
be separated. To proceed, it is crucial to decouple the rotation velocity and the wavelength 
to end up with expressions for the weak field approximation in terms of the observed fluxes. 
In order to do that, we assume that the shape of the local absorption profile is independent 
of the position on the stellar surface, apart from the Doppler shift introduced by v z , and 
that the CLV is wavelength independent. Furthermore, if we assume that the absorption 
profile is Gaussian, we end up with: 

1 /A-A o (l + ^(p,0)/c) 



/(A, v z (p, 9),fi) = /(/i)I(A, v z (p, 9),fi = l) = f(fji) 



r/exp 



2 V a 

(7) 

where I c stands for the continuum intensity, rj for the line depth, and a for the line width. 
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This approximation is not appropriate for stars with large radial temperature and density 
gradients because the lines close to the limb will differ strongly from the lines at the disk 
center. Apart from the Gaussian shape, it is possible to work with Voigt functions but the 
subsequent expressions turn out to be more complicated and no new physics is gained, so 
we prefer to stick with the Gaussian shape. The Gaussian shape is also consistent with the 
assumption of an unsaturated and weak spectral line that is necessary for Eqs. ($Z§ to hold. 

We can carry out the Taylor expansion of I(X,v z ,fi = 1) around v z — (we drop now 
the /^-dependence) to obtain the following expression for the intensity at each point of the 
stellar surface: 

r/A \ r/A „N 1 d n I(X,V 

n=l z 

Defining the emission line profile at zero velocity as 

•0(A) = r/exp 



v z =0 



1 ( A-_A 
2 



the absorption profile at an arbitrary velocity shift can be written as: 

1 d n I(X,v z ) 



I(X,v z ) = I c -^(X) + J2 



n=l 



(9) 



(10) 



v z =0 



Taking into account that the Gaussian function is the generating function of the 
Hermite polynomials: 



d n 
dx 1 



cxp 



x 

y 



-l) n exp 



x 



He n (x), 



where He n {x) is the Hermite polynomial of order n wi th weight function w(x) = ex pf— x 2 /2) 



also known as Chebyshev-Hermite polynomials (see 



Abramowitz fc Stegun 



19721 ). If a 



Voigt function is used inst ead for t he loc al profile, the derivatives are more complex and 



the recurrence relations of 



Heinzell ( 119781 ) can be used. Taking the previous definition into 



account, the derivatives (n > 0) of the intensity profile can be expressed as: 



d n I(X, v z 



dvl 



Vz=0 



ac I \ a 



(12) 



- 11 - 



so that the intensity profile of Eq. (J7|) is given by the following Gram-Charlier series: 



where we have extended the summation to n = to include the profile at zero velocity, 
since Heo(x) = 1. The previous formula shows that the intensity is given as the addition of 
infinite terms (potentially a large number of them for obtaining convergence) given by the 
product of a function that depends on the wavelength and another one that depends on 
the velocity along the line of sight. The latter is the only one that changes on the stellar 
surface, so that this variable separation allows us to carry out the surface integration. In 
order to simplify the following computations, it is interesting to define the function 



so that the intensity profile and its first and second derivatives with respect to the 
wavelength are given by: 



i{\,v„») = m[i c -^{\)u{\v z )\ 

I'(X,v z ,ii) = -^'(\)f( f i)'H(\,v z )-^X)f(fi)'H\X,v z ) 

/"(A, v z , fi) = -^"(A)/(/x)tt(A, v z ) - 2f(X)f(fx)H'(X, v z ) - ^(X)f(fi)n"(X, v z ), (15) 



where the ' and " symbols represent the first and second derivatives with respect to the 
wavelength. 

Applying the previous developments and dropping the dependences of the variables to 
avoid crowding, the polarized flux in the weak-field weak-line approximation can be written 



/(A, v z , /i) = /(/x) I c - ip(X) f^^ (-y He n (*—t*\ < 

n=0 n ' V°" C / \ ° J 



(13) 




(14) 
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as: 



F,(A) = F c -^fU 



C 



2 



^Q(A) = T A^G \f(\)(IP m - B*)fH + 2^(A)(S2 - B*)fH' + 1>{\)(B* - B%)fH" 
Fu(X) = ^-XtG [^"(X)B x B y fH + 2^{\)B x B y fW + ^{\)B x B y fU"\ 
F V (X) = C)fa[tf{\)BjH + 4(\)BjW\, (16) 

where we remind that the overlines represent integration in the plane of the sky, as indicated 
in Eq. (jl]). It is important to note that the area asymmetry of the polarized flux is zero as 
can be verified by computing the wavelength integral of the previous equations (it is easy 
to verify that the integral is zero by using the properties of the Hermite polynomials). 

The symbol F c stands for the continuum intensity flux and it is trivially given by: 

F c = 7Tl c (l - ~ - b -). (17) 



3 3, 

As expected, the Stokes parameters depend on quantities averaged over the stellar surface. 
Using the expansion on Hermite polynomials of the absorption spectral line, the terms 
appearing in the previous equations can be written as the following summations: 

n=0 N ' x ' 

n=0 v 7 v 7 

TO - 



n=0 

oo 



„ n\ \ ac / \ a 



13 



The functions Xi(n) are given by: 



x^ 


[n] = 


/ pdp 
Jo 


[ d0f(fi)B z (p,6 
'0 


K(P,6) 


X 2 \ 


[n) = 


/ pdp 
Jo 


I '* d6f(p)(B 2 x (p, 
lo 


d)-B 2 y (p,9)K(p, 




[n] = 


/ pdp 
Jo 


[ d9f(p)B x (p,6 
10 




X 4 I 


» = 


/ pdp 
Jo 


r ' d6f(pK( P ,6 

lo 


)• 



(19) 



It is important to note that Eqs. (fl6|) and (fT8|) are formally invariant for any magnetic 
field configuration, rotation profile and center-to-limb variation once the assumptions 
stated before (sufficiently weak magnetic field, wavelength independent CLV, unsaturated 
and weak Gaussian line) are fulfilled. The information about the magnetic configuration, 
rotation and the CLV is encoded in the Xi(n) functions, which may depend on additional 
external parameters. This is the case of the dipolar field with solid body rotation (see Sec. 
HJ), in which these functions depend on the orientation of the dipole and the angular velocity 
vector. 

The integrated Stokes parameters also depend on terms that include the wavelength 
derivative of the % function. Thanks to the effective separation between wavelength and 
velocity, these are easily obtained from the previous ones by using known properties of the 
Hermite polynomials. For instance: 

and the rest of terms are computed following the same strategy. 

When the star is not rotating, the only terms contributing to the summations in Eq. 
(TLB]) are those with n = 0. Likewise, the terms with derivatives of Eqs. fj20|) are always zero. 
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3. The weak-field approximation for observers 



The polarized fluxes of Eqs. (fT6l) can be used to explain stellar observations using a 
standard optimization method to fit the functional forms to the observations. However, 
we note that the expressions for the polarized fluxes depend on rj and a, properties of the 
local line profiles to which we do not have direct access. It is possible, though, to drop the 
dependence on r\ from the equations if we use the fact that ip'(X) and ip"{X) can be 

written in terms of Fj(X), F' T (\) and F'/(\), which are, oth erwise, observed quantities. This 



is the strategy followed by 



Martinez Gonzalez et al 



( 120121 ) to write the weak-field weak-line 



approximation for a non-rotating star in terms of the observed fluxes (and the parameters 
of the CLV and the Gaussian profile). So: 



(i/jfH' + ip'fn) 



F'/iX) = -(iPfU" + iP' f 7H + 2i,>fH>), 
which gives, after isolation of ip{X), ip'{X) and ^"(A): 

Fq(X) = 



(21) 



C 2 



XnG 



{Bl - B»)fH , J{Bl- Bl)fU> {Bl - BDfH fU> 



m 



^ F 'i 
m 



+ 



[bi - B$fw _ {Bi-Bpm fw (m - spm fW 



m 



2 {Bl-B*)U> fU> | 



m 



m 

■Fj) 



Fu(X) 



(■ 



B x B y fH , B x B y fW B x B y fH fW 
Fx + 2 



m 



==^— — Fj 



+ 



B x B v fH" B x B y fU fH" , BJ3jn fW B x B y H' fW , 

2 o 2 ±z^2 {r c - ti t 



m 



m 



Fy{X) 



-CXlg 



BJH BJHfW _BJW 



(22) 



These expressions are the generalization of the weak-field weak-line approximation to an 
unresolved star with an arbitrary velocity and magnetic field on the surface. The dependency 
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on the local absorption rj has disappeared once the observed Fj is used. However, the local 
broadening is still present and has to be either inferred from the observations or fixed to a 
reasonable value. The set of Eqs. (I22p represent a description of the Stokes parameters in 
terms of observed quantities and also in terms of wavelength-dependent unknowns that are 



explicit functions of the field and velocity distribution on the star (e.g., B z fH/fH). Both 
contributions are conveniently separated. In principle, one could use these expressions to 
obtain information about the magnetism of the star following two strategies. The first one 
is to parameterize these terms by fixing a model for the magnetic structure and infer the 
values of the parameters. This is the approach we follow in the next section. The second 
one is to non-parametrically infer these wavelength-dependent unknowns directly using a 
suitable inversion procedure. Since this would require additional regular izat ion techniques, 
a full treatment of the problem would have to be done in a Bayesian framework. 

An interesting property of the integrated Stokes parameters when rotation is taken 
into account is that the circular polarization flux depends on the derivative with wavelength 
of the flux (Fj) but also on the flux itself (Fi). Likewise, the linear polarization fluxes 
depend on the flux and its first and second order derivative with respect to the wavelength. 



When particularizing t he polarized flux of Eq. ([22]) to 



equations developed by 
equations. 



Martmez Gonzalez et al. 



;he non-rotating star, we recover the 



(120121 ) . This is a consistency check of our 



4. The rotating magnetic dipole 

Let us consider a star that is rotating with an angular velocity vector u), as shown in 
Fig. HJ In this section we use two different reference systems for our computations. The 
arbitrary reference system X'Y'Z, shown in the left panel of Fig. [2J is the one associated 
with an observer. The line-of-sight (LOS) is along the Z axis and the axis X' (arbitrarily 
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chosen by the observer) defines the reference direction for positive Stokes Q. The angular 
velocity vector is defined in this reference system by its modulus (u), the inclination angle 
(a) and the azimuth (7) with respect to the X' axis. The second reference system, XYZ, 
is obtained after a rotation of an angle 7 around the Z axis, so that the angular velocity 
vector is now contained in the XZ plane. This second reference system will help us in 
obtaining the final analytical expressions. 

One of the simplest non-trivial configuration of a stellar magnetic field is that of a 
dipole, whose axis is determined by a unit vector e along the dipole moment and an intrinsic 
dipolar field. The three components of the unit vector e on the cartesian system XYZ are: 

e x = sin i cos j 
e y = sin i sin j 

e z = cosz, (23) 

where i is the inclination of the dipolar moment with respect to the LOS and j is its 
azimuth with respect to the X axis. In the arbitrary reference system X'Y'Z, the azimuth 
is transformed to j' — j + 7. The emergent Stokes parameters Q(\) and U(X) are then 
transformed to the arbitrary r eference system follow ing Eq. fl3]). Note that our notation can 



be translated to that used by 



Landolfi et al. 



(119931 ) making the transformations (in format 



ours— )-theirs) : a — > i, j' — > A, 7 — y O, i — > I. Therefore, the cosine of the inclination of the 
dipole axis with respect to the rotational axis, usually termed (3, is plainly given by the dot 
product ijj ■ e, i.e., cos = sin a sin i cos j' + cos a cos i. 

The magnetic field at the surface of the star of a dipole placed at the center of the star 
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can be expressed in cartesian coordinates in the XYZ system as: 

3x(xe x + ye y + ze z )\ 
3y(xe x + ye y + ze z )\ 

3z(xe x + ye y + ze z )\ , (24) 

where B p is the intensity of the magnetic field at the poles. An arbitrary point (x, y, z) in 
the visible surface of the star can be written in terms of the polar coordinates on the plane 
of the sky £ as: 

x = p cos 9 
y = p sin 9 

Z = y/l - p\ (25) 

where we have selected the positive solution for z to isolate the visible half of the star. 

For simplicity, we assume that the star is rotating as a solid body because it greatly 
simplifies the equations (the introduction of a non-rigid rotation is formally identical, except 
for the computation of the integrals), so that: 

v = uj x r. (26) 

In terms of the polar coordinates on the stellar surface, the component of the velocity along 
the LOS in the XYZ reference system is given by: 

v z = up sin a sin 9. (27) 

Taking all these definitions into account and after some tedious algebra, it is possible 
to find relatively simple expressions for the X^n) integrals, which now depend on the set of 
parameters {i,j,a,cu,n,a,b}. The procedure is to plug the functional form of the dipolar 



B x 2 

B = -?*\e - 

B z = [e z - 



18 



i Z (LOS) 



* Z (LOS) 




-► 

Y 



Fig. 2. — This figure shows the two reference systems used in this work. We consider that the 
observer is along the Z axis, that is shared between both reference systems. The reference 
system X'Y'Z is an arbitrary one defined by the observer, where X' defines the reference 
direction for Q > 0. In order to simplify the computations, we use the XYZ reference system 
where the rotation angular velocity vector is in the XZ plane. The transformation between 
the two reference systems is given by a standard rotation with an angle 7. The vector e 
defines the dipolar magnetic field, while r defines the position on the stellar surface. 

magnetic field and carry out the angular and radial integrals on the visible disk. We end 
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These functions depend on the orientation of the dipole given by the angles i and j, 
the angular velocity parameters and the parameters of the CLV and the index n of the 
summation. The following functions take into account the dependence on the orientation of 
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the dipole and only depend on angular integrals: 

Ai(i,j,a,u),n) = S(n)cosi 
Bi(i,j,a,u,n) = —3S(n)cosi 
Ci(i,j,a,cu,n) = — 3S(n + 1) sini sin j 

A 2 (i,j,a,0J,n) = 3 cos 2 % [3S(n) — 8S(n + 2)] — 6 sin 2 % cos 2 jS(n) + 6S(n + 2) 
B 2 (i,j,a,uj,n) = 9 cos 2 i [S(n + 2) - S(n) + 2S(n + 4)] 

+ 9 sin 2 i cos 2 j [S'(n) - 4S(n + 2) + 4S(n + 4)] 

+ 9S(n + 2) - 18S(n + A) 
C 2 (i,j,a,cu,n) = 6 cos i sia i sia jS(n + 1) 
D 2 (i,j,a,uj,n) = 18 cos z sin z sin j [S(n + 1) — 2S(n + 3)] 
E 2 (i, j, a,uj,n) = sin 2 i cos(2j)S(n) 
A%(i, j, a,cu,n) = —3 sin 2 i cos j sin jS(n) 
B 3 (i,j,a,u,n) = 18 sin 2 i cos j sin j [S(n + 2) — S^n + 4)] 
C 3 (i, j, a, cu, n) = — 3cosisinicos jS{n + 1) 
D 3 (i,j,a,u,n) = 18 cos z sin i cos j [S(n + 1) — S(n + 3)] 
E 3 (i,j,a,u,n) = sin 2 2 cos j sin jS(n) 

A 4 (i,j,a,u,n) = S(n), (29) 
with the angular integral S'(n) given by: 

S(n) = jT d6 sin" 9=[1 + (-1)™] v^Fr T' 1 (l + |) , (30) 

which is valid for n > — 1. Likewise, the radial integrals -Fj(n) are independent on the 
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orientation of the dipole and are computed from: 
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(31) 



which are valid if n > — 2. 



Once the expressions for Xi(n) have been particularized to the case of a rotating dipole, 
it is possible to analyze in detail the shape of the profiles according to Eq. ( Tl6|) . We analyze 
the symmetry properties of the circular polarization profile, although this can be extended 
easily to the linear polarization profiles. Recall that the profile is a symmetric function 
of wavelength with respect to the rest wavelength of the line, while ^>'(A) is antisymmetric. 
According to the properties of the Hermite polynomials, the terms with n — 0, 2, 4, ... in the 



expression for B z fH in Eq. ([TBI are symmetric, while those with n — 1, 3, 5, . . . are always 



antisymmetric. Likewise, for B z fH', the terms with n = 0, 2, 4, . . . are antisymmetric, while 
those with n = 1,3,5, .. . are symmetric. Consequently, the polarized flux Fy(X) can be 
written as the addition of two functions, one symmetric and the other one antisymmetric: 

F V (X) = F^ m (X) + i^ ntisym (A), (32) 

where 



^ ym (A) 



j-iantisym / 



CX 2 g [*p'(X)B z fH antisym + ifj(X)B z fK' sym ] 
l (A) = CA^^(A)5jW sym + ^(A)5jW 7 a„t isym ]- (33) 
Taking into account the definition of the Xi (n) function of Eq. f|28|) , it is easy to verify that 
the symmetric contribution to -FV(A) depends only on sin i sin j, while the antisymmetric 



-22 - 



contribution depends only on cosi. Consequently, we will find standard antisymmetric 
circular polarization profiles whenever i — 0, while they will be fully symmetric when 
i = 90°. 

Stellar rotation is taken into account by keeping fixed the LOS and rotating the vector 
e around the angular rotation velocity uj. This is accomplished by computing the new 
rotated dipolar moment unit vector as: 

e rot = Re, (34) 

where the rotation matrix is given by: 

cos0 + sin 2 a(l — cos0) — cos a sin sinacosa(l — cos< 
cos a sin 
sin a cos a(l — cos 6) 



R 



- cos a sin < 

cos 
sin a sin d> 



cos • 



- sin a sin <p 

cos 2 afl — cos - 



(35) 



with indicating the stellar rotation phase and taking values between and 2ir. Note 
that our initial of phases is an arbitrary position of the vectors, while using the standard 
notation, the initial phase coincides with the uj and e vectors lying in the same plane 
parallel to the Z axis. This implies a phase shift between our <fi angle and the standard 
definition that is given by: 

sin i sin j 

tan O = . (36) 

cos a sin i cos j — sin a cos i 

Finally, the inclination and azimuth angles with respect to the XYZ reference system are 
obtained using the standard relations, taking into account that e rot is a unit vector: 

( e rot)j/ 



^rot = arccos [(e 



rot )z\ i 



jxot — arctan ■ 



e rot 



(37) 



5. Technicalities 



All the functions needed in the previous equations are well behaved and can be 
computed numerically without much complexity. The main difficulty resides on the fact 
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that, if the rotation velocity a; sin a is large, the convergence of the expansion in Eqs. (118]) 
is slow and oscillating, thus making it hard to compute. The main problem is that the 
Hermite polynomials grow in amplitude while the coefficients decrease slow. Only after 
many terms (the number of which depend on the ratio between the rotation velocity and 
the broadening), the coefficients start to compensate for the increase in the amplitude of the 
Hermite polynomials and convergence is eventually obtained. This is a well-known defect 



of the Gram-Charlier expansion (e.g. 



Blinnikov fc Moessner 



19981 ). When the rotation 



velocity is much larger than the Doppler width of the local spectral line (larger than a 
factor ~ 8), the compensation between the amplitude increase of the Hermite polynomials 
and the coefficients of the series expansion occurs for very large values of n [see Eqs. (ITS]) ]. 
At this moment, the double precision computations start to become unprecise and the 
terms in the expansion do not compensate adequately. This purely numerical drawback can 
be overcome if arbitrary precision methods are use d or if more elaborate metho ds like the 



Edgeworth expansion are adapted to our problem (IBlinnikov Sz Moessner 



19981 ) . In any 



case, we have been able to numerically compute expansions up to a; sin a ~ 8<jc/Ao. In 
other words, our method with a double precision summation of the expansions shown in 
Eqs. ({TBI) is presently applicable for rotation velocities not larger than 8 times the Doppler 
width of the local spectral line. 

In any case, it is advantageous to use the following set of modified Hermite polynomials: 

•AoV 1 



He* n {x) 



ac / n\ 



; He n (x) 



(3* 



which fulfill the following recursion formula: 



( -xHe n ^(x) 
ac I n 



A 



\ -L 



ac 



n 



He n - 2 (x) 



(39) 



Computing the Hermite polynomials with this recursion formula, the expansions of Eqs. 
(fT8"j) can be extended to much larger values of n. The advantage of using the modified 
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He^(x) Hermite polynomials is that their amplitude grow much slower with n than He n (x) 
due to the inclusion of the term Ao/(cxc) and n\. This facilitates the computation of each 
term in the expansions without numerical ove rflows. We note that the Her mite polynomials 



are computed using their recursion relations (lAbramowitz fc Stegun 



1972). 



From the computational point of view, once the rotation angular velocity of the star is 
set, the wavelength dependence of Eqs. (jl8J) and their derivatives in terms of the Hermite 
polynomials can be precomputed once and used in all subsequent calculations. This greatly 
accelerates the process because only the Xi(i,j,n) functions have to be computed for a 
given orientation of the dipolar moment and added until convergence. 

A free computer program in IDL is available to compute the observed Stokes fluxes 
for a rotating dipole using the expressions in this papej^. In this non-optimized version of 
the code, we synthesize profiles evaluated at 500 wavelength points with a step of 8 mA in 
~35-40 ms. 



6. Illustrative examples 



Given our assumptions (weak-field, weak-line, fixed center-to-limb variation, global 
dipolar field), one of the most direct applications of our equations is to study the magnetism 
of hot stars. For t his reason, we use ty pical values of the broadening and rotation velocity 



of such stars (e.g. 



Petit fc Wade 



20121 ). We show in Fig. [3] two examples for a dipolar 
magnetic field with different field configurations. Both computations share the parameters 
for the CLV (a = 0.2 and b = 0.2), the magnetic field strength in the pole (B p = 500 G) and 
the rotation velocity on the surface (v = 40 km s _1 ). The artificial selected spectral line is 
centered at Aq = 5000 A and has an effective Lande factor of g — 3 and its equivalent for 



2 http : //www. iac . es/project/magnetism/dipole 
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linear polarization of G = 9. Since the magnetic field takes its maximum value at the poles, 
the spectral line can be safely considered in the weak-field regime. The line broadening is 
set to 0.1 A, typical of an observation at a spectral resolution of R = 65000 of a line with a 
microturbulent width of ~5 km s _1 . These values of the broadening and rotation velocity 
are close to the numerical limit that we presently have for computing the summations of 
Eqs. ( I18p . Given this limitation, our scheme is presently not applicable to fast rotators 
with narrow local absorption profiles. The computed intensity flux normalized to the 
continuum flux is shown in the upper panel of each figure. It is independent of the phase 
since, under the weak field regime, the intensity does not depend on the magnetic field and 
we have assumed that the star is structureless. However, it takes into account the rotation 
broadening and the presence of a CLV. Concerning the linear polarization fluxes, we see 
maximum amplitudes around 0.02%, while this increases to values around 100 times larger 
for circular polarization. This is a direct consequence of the fact that linear polarization 
appears only in second order. 



7. Conclusions 

We have generalized the expressions for the Stokes parameters in the weak-field 
weak-line regime to the stellar case in which rotation is taken into account. This produces 
an intensity flux profile that is broadened by rotation while it is insensitive to the magnetic 
field. As a consequence, our formalism can only be applied to stars in which the intensity 
profile is not observed to be modulated by rotation. However, the circular and linear 
polarized fluxes are sensitive to the topology of the magnetic field. We have demonstrated 
that the inclusion of Doppler shifts generates a dependence of circular polarization of 
the intensity flux itself, and not only on its wavelength derivative. Likewise, the linear 
polarization profiles depend now on the first and second wavelength derivative and on the 
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intensity flux itself. We have verified also that the expression s developed in this paper 
simplify to those obtained by 



Martinez Gonzalez et al 



(J2012J) in the non-rotating dipole 



case. The general weak-field expressions have been particularized to a rotating dipole. The 
extension to a more general multipolar magnetic field is relatively easy and can be done 
fast if analytical expressions for the angular and radial integrals exist. In subsequent papers 
of this series we will deal with thin magnetized disks in Keplerian rotation and unresolved 
structures in the quiet Sun. Likewise, we will also explore the possibility of using more 
elaborate expansions to improve the convergence. 



We thank R. Manso Sainz for fruitful discussions. Financial support by the 
Spanish Ministry of Economy and Competitiveness through projects AYA2010-18029 
(Solar Magnetism and Astrophysical Spectropolarimetry) and Consolider-Ingenio 2010 
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Wavelength [A] Wavelength [A] 



Fig. 3. — Examples of the flux normalized to the continuum intensity flux for a full rotational 
phase. Both examples share a = 0.2, b = 0.2, B p = 500 G and u>R star = 40 km s" 1 , where 
-Rstar is the stellar radius. The left panel corresponds to the case % = 45°, j' = 80°, a = 70° 
and 7 = 70°. The right panel has i = 45°, j' = 30°, a = 90° and 7 = 20°. All polarized 
fluxes are given in percentages. 
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